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Abstract. The most general vortex solution of the Liouville equation (which arises in non- 
relativistic Chern-Simons theory) is associated with rational functions, f(z) = P(z)/Q(z) 
where P(z) and Q(z) are both polynomials, degP < degQ = N. This allows us to prove 
that the solution depends on 4N parameters without the use of an index theorem, as well 
as the flux quantization : $ = — 4A%(sign «). 
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1. Vortex solutions 

Non-relativistic Chern-Simons theory supports vortices [1]. (See also [2] for reviews). 
For suitable values of the parameters, these vortices arise as solutions of the first-order 
"self-duality" (SD) equations, 

(1.1) (£>i ±iD 2 )V> = 0, kB = -q, 

where ip is a complex scalar field minimally coupled to a static gauge potential A = (A x , A x ) 
in the plane; Djip = (dj — iAj^tp is the covariant derivative, B = d x A y — d y A x is the 
magnetic field, and g = ip*ip is the particle density. The solutions we are interested in are 
"non-topological" in that g — > at infinity. 

Using the complex notations d = \{d x — id y ), A = A x — iA y , z = x + iy, the first SD 
equation are solved by A = —2idlnip. Setting ip = g 1 / 2 e luJ , we get 

(1.2) A = -idlng + 2du. 

Reinserting A into the second SD equation, we end up with the Liouville equation, 

2 

(1.3) Ahi£=- — g, 

\k\ 

whose all real solutions are known. They are in fact given in terms of an arbitrary analytic 
function on the complex plane, 

\f'(z)\ 2 

(1.4) e = 4|«' ' ' 



(1 + I/WI 2 ) 2 ' 

For f(z) = z~ N , for example, we get the radially symmetric solution 

r -2(N+l) 



(1.5) q = 4\k\N' 



(l + r~ 2N ) 2 



To be regular, has to be here an integer at least 1 [1], [2]. 

A striking feature of these non-topological vortex solutions is that their magnetix flux 
is an even multiple of the elementary flux quantum, 



(i.6) *=y 



$ = J Bd 2 x = 2N x $ , $ = -2tt, AT = 0,±1,... 

This can be seen easily in the radial case, taking into account the asymptotic behaviour, 
g oc r - 2 ( N + 1 ) as r — > oo, of the particle density. We are not aware of a general proof, 
though. 
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Another peculiarity of these vortices is that, for fixed N, the solution depends on 4N 
parameters. The proof given by Kim et al. [3] uses the Atiyah-Singer index theorem. If 
such an approach is perfectly justified for "Nielsen- Olesen" vortices or "BPS" monopoles 
[4], it seems too-powerful when, as in our case, all solutions are known explicitly. 

A 4A r -parameter family of solutions can be written down at once: consider 



Plotting the particle density allows us to interpret the associated solution as representing 
N separated vortices located at the points z i: with individual scales and phases c; [5]. 
With some work, the associated flux is found to be $ = — 4A r 7r(sign k) [1]. 

The questions we address ourselves sounds: what is the most general / yielding phys- 
ically admissible solutions ? What is the associated magnetic flux ? Can we count the 
parameters without an index theorem ? Below, we show that imposing suitable regularity 
conditions only allows rational functions. In detail, we prove 

Theorem 1 . Let us consider a vortex solution of the Liouville equation with finite magnetic 
flux $ < oo. If the magnetic field B is regular on the complex plane and is such that r 2+s B 
(S > 0) is bounded when r — > oo, then f is a rational function, 



where P(z) and Q(z) are polynomials with degP < degQ. The coefficient of the highest- 
order term in Q(z) can be normalized to 1. 

It follows that the general solution does indeed depend on 4N parameters, namely on 
the 2 x N complex coefficients of the polynomials with P(z) and Q(z). The particular 
form (1.7) is recovered by expanding / into partial fractions, provided Q(z) only has simple 
zeroes. 

Next we prove 

Theorem 2 . The magnetic flux of the vortex associated with the rational function f(z) in 
Eq. (1.8) is $ = — 4A r 7r(sign k), where N is the degree of the denominator, 



N 



(1.7) 





(1.9) 



N = degQ(z) 



and is hence an integer. 

As a corollary, we also get the general theorem on flux quantization, as in Eq. (1.6). 
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We establish our theorems by elementary complex analysis [6], [7], as the result of a 
series of Lemmas. We start with proving the 

Lemma 1 : Let f(z) be a complex function which only has isolated singularities. Let 7 be a 
curve in the complex plane which avoids the singularities of f ; set zo = f(0) and zi = f(l). 
Then 

!/(*,) -/(*,)! < / 
Vi + i/(^)lVi + l/(--i)l 2 " Ai + l/WI 2 

where \dz\ = \d"f/dtdt\. 

This proposition has a nice geometric meaning: the left-hand side is the length of 7 
with respect to a metric inherited by stereographic projection, while the right-hand-side is 
the integral of g 1 / 2 , the square-root of the particle density. 

Next, using Lemma 1, we demonstrate 

Lemma 2 . The function f can not have an essential singularity in the complex plane. 

This lemma eliminates functions like f(z) = e 1 / 2 . This function would yield in fact a 
density which is unbounded at the origin. Writing z = r exp (i9) , the density becomes 

4exp(^) 



r4(l + exp (^)) 



2 ' 



so that on the imaginary axis, 9 = tt/2, lim g = +00. 

i — >o 

Next, we prove 

Lemma 3 . The function f can not have an essential singularity at infinity. 

Lemma 3. rules out the functions f(z) as e z . For this choice, the particle density 
reads 

4 

e = 72- 

exp (— r cos 9) + exp (r cos 9) 



On the imaginary axis 9 = tt/2, we have g = 4: the particle density is not localized. 

Now a theorem found in Whittaker and Watson [6], 5.64., p. 105. says that the only 
one-valued functions which have no singularities, except poles, at any point (including 00), 
are rational functions. This allows us to conclude that our f(z) is indeed rational. Then, 

Lemma 4 . The polynomials in Eqn. (1.8) can be chosen so that degP < degQ, the highest 
term in Q having coefficient equal to 1. 

This proves our Theorem 1. 
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Now, to evaluate the magnetic flux of the vortex associated with the rational function 
/ in Eq. (1.8), we show that 

Lemma 5 . Let z\, . . -Zn q denote the distinct roots of the denominator Q(z), each having 
multiplicity n^. The particle density (1-4) can also be written as 



(1.11) 



q = 4\k\8 



( df, |/| 2 
1 / } l + l/l 2 



Nq 

+E 

i=i 



Hi 



Z - Zi 



where the bracketed quantity is a regular function on the plane. 

Then the flux (1.6) is converted into a contour integral at infinity, by Stokes' theorem, 



(1.12) $ = 2z(sign k) j> 



■df, |/| 2 
- / ] 1 + l/P 



Nq 



^ Z-Zi 



i=l 



dz, 



where S = Soo is the cercle at infinity. 

The integrand of (1.12) is related to the vector potential. Using (1.2), this latter reads 
in fact 



(1.13) 



. .,df, |/| 2 - ( 9 2 f^ 9 



/ y i + l/l 

To get a regular A, the phase u has to be chosen so that 



rii — 1 \ - m,- — 1 



, z - z l 
i=i t=i 



so that the integrand in (1.12) is 



Nq N P 

rii + 1 \ -v mi — 1 



M^)+<E^-E^ 



i=i 



i=i 



the integral of the terms in the curly bracket on the circle at infinity vanishes. 
Now using 



Lemma 6. 



(1.15) 



3 [ f'l + l/l 2 ' 



6 



Horvdthy & Yera 



the second term in (1.12) is evaluated at once, 

„ n q n q 
(1.16) $ = 2i <b y^-^dz = -4tt(YV) = -4iV7r(sign k), 

J b 1=1 % i=l 

where N is the degree of Q(z). This yields Theorem 2. 

The flux has been previously related to the inversions [8]. Their argument goes as 
follows: the particle density behaves as p ~ r z(N-i) w h e n r — > 0. The regularity of 
the vector potential requires the phase to be chosen as uj = (N — 1)9. Then the inver- 
sion symmetry implies the behaviour p ~ r -2(JV+i) at infinity, so that the flux is indeed 
2n(N - 1) + 2yr(iV + 1) = 4nN. 

This argument is only valid in the radial case, though. To see this, let observe that 
the choice f(z) = (1 + z)~ 2 - 2(z - l) -1 , yields, for example, flux $ = —12n (i.e. = 3) 
and is interpreted as a 2- vortex sitting at z = — 1 and a 1- vortex sitting at z = 1. The 
particle density does not behave as claimed by Kim et al [8], rather as p ~ r 6 (instead of 
r 4 ) when r — > 0, and as p ~ r -4 (instead of r -8 ) when r — > oo. 

Where does the error come from ? On the one hand, the behaviour at the origin 
assumed by Kim et al. is consistent with our formulas (1.4), (1.14) in the radial case 
only. On the other hand, albeit the Liouville equation is indeed inversion-invariant (indeed 
invariant with respect to any conformal transformation) this is not true for individual, 
non-radial solutions. Therefore, the large-r behaviour of a solution can not be inferred 
from that for small r. 

In this paper, we only considered the case of single- valued functions /. It seems 
however, that multiple-valued function do not qualify. For example, the charge density 
g associated with f(z) = In z, is also multiple-valued, and hence physically inadmissible; 
remember also that in the radial case f(z) = z~ N , the regularity requires iV to be an 
integer. 
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2. Proofs. 



Proof of Lemma 1 . Stereographic projection carries over the natural metric from the Rie- 
mann sphere to the complex plane. The scalar product of two tangent vectors, u ad v, at 
a point p of the plane is 



(2.1) 



9 P (u,v) 



(i + M 2 ) 2 



u ■ V, 



where "•" is the ordinary scalar product in R , \p\ 2 = p ■ p. The length of a curve 
[0, 1] 3 t h-> F(t) G C w. r. to this metric is 



(2.2) 



L(T) 



f 


oT 




dT 


Jo 


dt 


dt, 


dt 



dr dr\iV 2 



l + \T(t)\f 



Then the distance of two points, wq et wi, in the complex plane is the l'infimum of length 
of the curves between the points, 

d(w ,w 1 )=Tn£ r {L(T) | r(0) = w , T(l) = Wl }. 

Let us now consider an analytic function w = f(z). f can also be viewed as a mapping 
of the z-plane into the w-plane; the latter is endowed with the distance defined here above. 
If 7(t) is an arbitrary curve in the z-plane with end-points zq and z±, its image by / is a 
curve T = fo"f in the w-plane with end-points wq = f(zo) and w\ = f(zi). By (2.2), the 
length of r is the r. h. s. of (1.10), 



(2.3) 

Thus 
(2.4) 



2|/'(s)| 



+ \f(z)\ 2 



\dz\. 



d(w ,w 1 ) = d(f(z )J(z 1 )) < 



1 + 



\dz\ 



But the distance on the w-plane is just the distance on the Riemann sphere. But this 
latter sits in R 3 , so that the (geodesic) distance on the sphere is greater or equal to the 
natural distance in R 3 : 



(2.5) 



d(w ,w 1 ) > 







- Wi 






Wq 


2 V 1 + 




2 



equality being only achieved for w± = wq. Setting Wi = f(zi), % = 0, 1, the inequality 
(1.10) is obtained. 

Q. E. D. 
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Proof of Lemma 2 . Let us assume that / has an isolated essential singularity at a point zq. 
Then it is analytic in some disk D = D(zq; e) \ {zq}. Now, according to Picard's Theorem 
([7], p. 90): If zq is an isolated singularity of a holomorphic function f(z), then for each 
r > 0, I'image of the annular region {z e C|o < | z — zq\ < r} is either the whole of C or 
C without a single point. 

Let us first assume that z\ is a point in D such that f(zi) = 0. Then, since the 
particle density, g, is a regular function on the plane which goes to zero at infinity, there 
is a real number M such that 



4|/'(*) 
(l + l/WI 2 ) 1 



= Tr^rrkv2 < M *> Vz e C. 



Eqn. (1.10) in Lemma 1. with zq = z, yields, using f{z\) = 0, 



<M I \dz\ 

J <y 



for all curve 7 s. t. 7(0) = z et 7(1) = z\. For the straight line 7(t) = z + t{z\ — z) in 
particular, the r. h. s. becomes \z\ — z\M < 2eM, since z\ and z both belong to D. Thus, 
chosing e to have 4Me < 1, \f(z)\/y/l + \f(z)\ 2 < 1/2, which implies that \f(z)\ < 1. The 
function /(z) is hence bounded in D, which contradicts hypothesis that zq is an essential 
singularity. 

Now if / does not vanish in D = D(^o;e), one can chose z' with /(^o) sufficiently 
small so that f(z) is still bounded in D. 

Q. E. D. 

Proof of Lemma 3 . Let us assume, on the contrary, that infinity is essential singularity of 
f(z). From the large-r behaviour B = o(r 2 ) we shall deduce that / is again bounded at 
infinity, a contradiction. 

In fact, if 00 is an isolated singularity, then f(z) is holomorphic in some neighbourhood 
D = {z G C \z\ > iV}, of infinity. In this neighbourhood, one can find a point Zq where 
/ vanishes f(zo) = by Picard's theorem ( 1 ). 

Now, due to the imposed groth condition on B = —q/k, N can be taken so that 

(2.6) q<C 2 \z\- 2 ~ s , for all \z\ > N, 

where C > is a constant. Let us chose iV such that 4C(tv + -) < N 6 ^ 2 . Then, for all 



complex number z, \z\ > \zq\, f is bounded, \f(z)\ < 1. To see this, consider zi, the 



( x ) If f{z) never vanishes in the region D, it is enough to chose f(zo) small enough. 
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intersection of C, the cercle around with radius \z\, with the straight half-line [0, zq). 
Then 

(2.7) d(0, /(*)) < d(0, /(zi)) + d(/(zi), /(*)). 
Now, applying the inequality (2.4) to the circular arc (zi, z) we get 

(2.8) d(/(*i), /(*))< [ e {f{z)) 1/2 \dz\. 

J arc 

From this we deduce, using (2.6), that 

(2-9) d(f(z 1 )J(z))< WJ ^. 

On the other hand, applying (2.4) to the segment "f(t) = zq + t{z\ — zo), t e [0, 1], we get 

i/2, , , f C . , . 2C 
I 1 \dz\ < / 1 " 



(2.10) d(0,/(zi)) < | ^(/(z)) 1/2 |d^l < | uj^l^l < 



2 



when the condition (2.6) is used again. The inequalities (2.7)-(2.9)-(2.10) imply that 

2(1 1 1 

(2-11) d(0,/(z))<^(7r+^)<-; 

Now, by (2.5), we have 

d(0,f(z))> lf{z)l • 



Vl + l/WI 2 ' 

using (2.11), we get finally f(z)\ < 1. The function f(z) is hence bounded in some 
neighbourhood of infinity, so that this point cannot be an essential singularity. 

Q. E. D. 

Now (as explained in Chapter 1), Theorem 5.64 of Whittaker and Watson [6] allows 
us to deduce that / is a rational function, f(z) = P(z)/Q(z), where P(z) and Q(z) are 
both polynomials. 

Proof of Lemma 4 . Now, since / and / _1 are readily seen to yield the same solutions, 
we can assume that degP < degQ. The case degP = degQ is eliminated by a simple 
redefinition, as in Ref. [1]. In fact, 

m = fo + B^y 

where A(z) are B(z) polynomials s. t. degA(z) < deg B(z) and fo ^ 0, is readily seen to 
yield the same density (1.4) as 



B(z) 
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with the polynomials A(z) and B{z) defined as 

Then the coefficient of the highest-order term in Q(z) can be normalized to unity. 



Q. E. D. 



Proof of Lemma 5 . In complex notations, the general solution (1.4) is expressed as 



(2.12) 



q = 4\k\8 



df, \f\' 



where 3 = \ [d x + idy) . Let us now denote by Zi, i = 1, . . . , Nq the distinct roots of the 
denominator, Q(z), each having a multiplicity rii. [Since P(z) and Q(z) have no commun 
roots, these are the same as the the poles of f(z)]. Then the function 



N C 



(2.13) 



Hi 



( df, |/| 2 ^ 

1 / j i + m 2 + ^-^ 



is regular on the complex plane. Indeed, in the neighbourhood of a root, z^, of order rii of 

Q(z), 



f 



(z-ZiY 



df l/l 2 
/ ~ (z-ZiY 1 + |/P 



1. 



In contradistinction, in the neighbourhood of a zero (denoted by Zq) of order k > 1, 
of P(z) [which is the same as a zero of f(z)], we have : 



/ ~ (z - Z ) k 



df k 



\f\'< 



~ \ z — z o\ j 



/ z-Z ' 1 + |/| 2 



so that (2.13) is regular. Now, since 



the second term in (2.13) can be added to the expression (2.12) of the density, yielding 
(1.11). 

Q. E. D. 
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The integrand in (2.13) being a regular function, the flux $ = —(1/k) J gd 2 x can be 
converted, by Stokes theorem, into an integral at infinity, to yield (1.12). 

Proof of Lemma 6 . Let us now denote by Zi, i = 1, . . . , Np, the distinct roots of the 
numerator, P(z) [and hence those of f(z)]. Using the theorem on the residues, 

since lim f(z) = 0. 

Q. E. D. 
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